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Summary. We sketch out a new geometric framework to construct Hamiltonian operators 
for generic, non-evolutionary partial differential equations. Examples on how the formalism 
works are provided for the KdV equation, Camassa-Holm equation, and Kupershmidt's defor- 
mation of a bi-Hamiltonian system. 



1 Introduction 

In this short paper we will discuss the following question: What happens to a Hamil- 
tonian operator of an evolution system if we change coordinates so that the system 
becomes non-evolution? 

Using the traditional definition of a Hamiltonian structure one cannot answer this 
question, since the definition is tied to evolution form of the system at hand. However, 
first, not all equations have a natural evolution form, and, second, an evolution form 
of a system of equations is not unique. Let us consider some examples. 

Example 1 (KdV). It is well known that the KdV equation — u^xx + (iuu^ has two 
compatible Hamiltonian operators: 

so that the equation can be written in the following ways: 

Mf = Uxxx + 6mMj = Dx^{u - 
ou 
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= {Dxxx + ^uD, + 2u^) ^ (m72) , 
ou 

where 8/5u denotes the Euler operator (the variational derivative) and is apphed to 
the two Hamihonian densities. 

Let us introduce new dependent variables v and w and rewrite the KdV equation 
in the form 

Ux = V, Vx = w, Wx = Ut — 6uv. (2) 

In the new coordinates, the KdV still has an evolutionary form, but with respect to 
another independent variable (x instead of t). A natural question arises then: Is the 
KdV equation in the form (2) Hamiltonian? An affirmative answer to this question 
was obtained by Tsarev in [9]. He proved that transformations of the type (2) preserve 
the Hamiltonian property of all evolution systems for which the Cauchy problem is 
solvable. Our approach is very different from Tsarev's one. Below we explain why 
this fact holds true for all transformations of variables and without the assumption on 
the Cauchy problem. We will also show how to compute the Hamiltonian structure 
in new coordinates. For the above example the answer is the following: 



u 



-1 








—6m 


6m 









2m 






{uw-v^/2 + 2u^) 



-2m -Dt-2v \ /5/5m^ 

A -I2u^-2w\ 5/5v I (-3mV2-w/2)- 



Dt + 2v I2u^ + 2w SuDt + 4ut J \5/5wj 



(3) 



Example 2 (Camassa-Holm equation). Camassa and Holm have written their equa- 
tion Ui — Utxx — uUxxx — '^UxUxx + 3uUx = in a bi-Hamiltonian form by introducing 
the new variable m = u — Uxx- The equation now takes the form 



with 



nit = -uiUx - 2uxm = Bi — — = B2 — 
dm dm 



Bi = —{mDx+Dxm), M'\ = \ \ mudx, 



(4) 



2 

B2=dI-Dx, Jf2 = ^J{u^ + uu^}dx. 

Note that ^ and are viewed as functionals in m and u, but not in m solely. To 
get rid of m, one is forced to assume that m = (1 —Dx)^^m in the Hamiltonian densi- 
ties. The use of the inverse of the operator 1 — is not elegant from mathematical 
viewpoint. We will find a bi-Hamiltonian structure for the Camassa-Holm equation 
written in the initial non-evolution form and thus get rid of the term (1 — 
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Example 3 (Kupershmidt deformation). Consider a bi-Hamiltonian evolution system 
of equations m, = f{t.x, u,u^, . . .), u and / being vector functions, with compat- 
ible Hamiltonian operators A i and A2 and a Magri hierarchy of conserved densities 
Hi, H2, . . . 

Dt{Hi)=0, Ai— — =A2^ — . 

OU OU 

In [8J, Kupershmidt defined what he called the nonholonomic deformation of the 
above system: 

M,=/-Ai(w), A2(w)=0. (5) 

We call system (5) the Kupershmidt deformation of the system ut = f. The motivating 
example of this construction is the so-called KdV6 equation (see [4J) 

Ut = Uxxx + 6uUx-Wx, Wxxx + ^UWx + 2UxW = (6) 

which is the Kupershmidt deformation of the KdV equation. The authors of [4] have 
shown that the KdV6 passes the Painleve test and conjectured that the system is 
integrable. Kupershmidt, in [8], found a hierarchy of conservation laws of the KdV6 
as a particular case of the following general fact. 

Theorem (Kupershmidt). Let Ut = f he an evolution bi-Hamiltonian system, with 
Ai, A2 being the corresponding Hamiltonian operators. If this equation has a Magri 
hierarchy of conserved densities ^ = 0, Aj ^ = A2 then Hi, H2, . . . are 
conserved densities for (5). 

Proof. 

Kupershmidt also conjectured that Hi, H2, ... commute in some sense so that the 
KdV6 is indeed integrable. Below we will see that this is true and, moreover, sys- 
tem (5) is bi-Hamiltonian. 

Our framework to study Hamiltonian structures for general PDEs is the geome- 
try of jet spaces and differential equations. We assume the reader to be familiar with 
the geometric approach to differential equations and hence we include only the no- 
tation and the coordinate descriptions in the next section. We refer the reader to the 
books [1,6] for further information. 

2 Notation: infinite jets and differential equations 

In what follows everything is supposed to be smooth. 
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We denote an infinite jet space by J°°. This can be the space of jets of submani- 
folds, maps, sections of a bundle, and so on, and it is not important to us here. Coor- 
dinates on y~ are x; (independent variables, i = 1 , . . . ,n) and (dependent variables, 
j = I,. . . ,m, a being multi-indices). 

The formulas 

^ d ^ j d 



provide expressions for the total derivatives. The vector fields D, span the Cartan 
distribution on 7°°. To every vector function on 7°°, there corresponds the evolutionary 
field 

j,a oUa 



The matrix differential operator 



ydf_ 



is the Unearization of a vector function /. It is defined by the formula ^/((jo) = E(p{f). 
The linearization is a differential operator in total derivatives; we shall call such 
operators ^-differential operators. 

The coordinate expression for the adjoint ^-differential operator is 



A* = 



if A 



Let Fi^{xi,Uff) = 0, k = 1, . . . /, be a system of differential equations. Then the 
relations F — {Fi,...,Fi) = together with Dff{F) = define its infinite prolonga- 
tion # c For the sake of brevity we shall call the infinite prolongation of a system 
of differential equations the equation. The operator £s = £f\s is the Unearization of 
the equation S". 

In this paper, we only consider equations # whose linearization £^ is normal in 
the following sense. 

Definition 1. A ^-differential operator V called normal if the compatibility opera- 
tors for both V and V* are trivial. In other words, if there exists a ^€ -differential 
operator A such that A oV = on S' then A =0 on S' as well, and the same holds 
true with V* instead ofV. 

An evolutionary field Eq, is a symmetry of the equation if E(p{F)\g = £s{<p) = 
0. If E(p is a symmetry then (p is said to be its generating function. We often identify 
symmetries with their generating functions. 

A vector function S = {S^,...,S") on S' is a conserved current if L/^K'^') = 
on (o. A conserved current is trivial if there exist functions 7]^ on such that S' = 
Lj<i 
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Conservation laws of S are classes of conserved currents modulo trivial ones. 
To every conservation law, there correspond its generating function, which is com- 
puted in the following way. If S = (5 5") is a conserved current, so that 
1^, D,(5') = on then there exists a ^-differential operator 4 such that J^, A (5') = 
A (F) on J°°. The generating function of the conservation law is defined by = 
{xj/i,..., Xj/m) =A*{1). Note that i/a = if and only if the conserved current S is triv- 
ial. One can prove that every generating function y/ satisfies the equation (i/a) = 0, 
so that the set CL{S') of conservation laws of ^ is a subset in the kernel of i*^, 
CL(#) c kerl*^. 

3 Cotangent bundle to an equation 

Let us introduce our main hero. For every differential equation S' we define a canoni- 
cal covering T* : ^*(#) called the t -covering. The equation ^*(#) is given 
by the system 

tp{p)=Q, F = 0, 

if is given by F = 0. Here p = {p^ , . . . are new dependent variables, I being 
the number of equations F = (Fj, . . . We endow with the structure of a 

supermanifold by choosing the variables p^ to be odd. The covering t* is the natural 
projection T* : (m^^Pct) ^ ("cr)- 
Note that 

{F,p)^t,Fip' (7) 

i=l 

is the Lagrangian for the equation 

It is easily shown that £ ^*(g) is normal if £g is normal. 

From the above definition it is not seen why we said that £*-covering is canonical. 
Indeed, the definition uses the embedding S" ^J°°, but later we will show that ^* ((#) 
is independent of the choice of this embedding. 

Remark 1. For an arbitrary -differential operator A one can define the /i-covering 
in the same way as the £* -covering is associated with the operator 

The most interesting for us property of the ^*-covering is given by the following 
theorem. 

Theorem 1. There is a natural 1-1 correspondence between the symmetries of S and 
the conservation laws of^*{S') linear along the fibers ofT*. 

The expression "linear conservation law" means that the corresponding conserved 
current is linear along the fibers of T* (i.e., linear in variables p'^). Here and below 
we skip the proofs that can be found in our joint paper with S. Igonin [3]. Let us 
nevertheless describe the correspondence stated in the theorem in terms of generat- 
ing functions. If (jj is a symmetry of equation S' then there exists a "^-differential 
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operator A such that £f{(p) = A{F). Consider the adj oint operator A*. It can be natu- 
rally identified with a fiberwise linear vector function (jO/j on^*{S'). Then the vector 
function (9, (jO/j) is the generating function of the conservation law that corresponds 
to the symmetry (p. 

In the geometry of differential equation it is very useful to construct an analogy 
with geometry of finite dimensional manifolds. We shall now use this approach to 
clarify the meaning of the above theorem. Let us start building our analogy with the 
following two rather standard correspondences (cf. [10] and references therein): 



Manifold M 



functions 
vector fields 



conservation laws 
symmetries 



Now, using Theorem 1, we can say that the analog of the -covering is a vector 
bundle such that vector fields on the base are in 1-1 correspondence with fiberwise 
hnear functions on the total space of the bundle. Obviously, such a bundle is the 
cotangent bundle. So, the f *-covering is the cotangent bundle to an equation, and we 
can continue our manifold-equation dictionary: 



Manifold M 

functions 
vector fields 
T*{M) 



conservation laws 
symmetries 



Remark 2. This dictionary can be easily extended: 



Manifold M 

functions 
vector fields 

T*{M) 
T{M) 

De Rham complex 



PDE<r 

conservation laws 
synometries 



,l,n-l 



^2,«-l 



Here ^{S) is the ^-covering (see Remark 1). The complex E-^ 



,0,n-l 



' • is (« — l)st line of the Vinogradov '^^'-spectral sequence (see [10] and 
references therein). In this paper we use only the first three entries of the dictionary. 

Remarks. In [7], Kupershmidt defined the cotangent bundle to a bundle. This con- 
struction can be identified with the ^*-covering of the system 



u] = 0, 



uj = 0, 



0. 
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At this point, a natural question may arise: what is the analog of the Poisson 
bracket on the cotangent bundle? The answer is that the ^*-covering is endowed with 
a canonical Poisson bracket. More precisely, since we changed the parity of fibers 
in the £*-covering, this bracket is a superbracket and is the analog of the Schouten 
bracket. We shall call it the variational Schouten bracket. 

To define the bracket, recall that * ) has the Lagrangian structure (7). Hence, 
by the Noether theorem there is a 1-1 correspondence between conservation laws 
on ^*((?) and Noether symmetries of ^5f*(f5'). If xj/ is the generating function of 
a conservation law, then is,^ is the corresponding Noether symmetry. The set of 
Noether symmetries is a Lie superalgebra with respect to the conamutator, so we 
obtain a structure of Lie superalgebra on conservation laws on ^*((#) uniquely de- 
termined by the equality 

According to our manifold-equation dictionary, conservation laws on ^* (^) cor- 
respond to functions on T*{M). The latter are skew multivectors on M (this is why 
we have changed the parity of fibers of the £*-covering — to get skew- symmetric mul- 
tivectors). So, we shall call conservation laws on ^*{S') the variational multivectors. 
Linear conservation laws, as we saw, are vectors, biliner ones are bivectors and so on. 

The generating function of a variational /c- vector is a vector function on 
which is {k— l)-linear along T*-fibers. Such a function can be identified with a 
(k — 1 ) -linear '^'-differential operator on S . In coordinates, this correspondence boils 
down to the change pa ^ D^. Thus, we can (and will) identify variational multivec- 
tors to multilinear "^-differential operators. 

More precisely, in the above identification we will use not operators but equiv- 
alence classes of "^-differential operators modulo operators divisible by £*g. This 
is being done, because operators of the form Ooi*^ correspond to trivial functions 
on ^* (ff ) . But we wiU not change terminology, we say operator instead of the equiv- 
alence class. 

For the sake of brevity and because we are interested in the Hamiltonian for- 
malism, let us restrict ourselves to bivectors, which are identified with linear 
differential operators. Formulas presented below for bivectors (= linear operators) 
can be easily generalised to multivectors (= multilinear operators). 

Theorem 2. An operator A is a variational bivector on equation § if and only if it 
satisfies the condition 

lsA=A*tg. 

Remark 4. If S' is written in evolution form then the above condition impUes that 

A* = -A. 

From this theorem it follows that a Hamiltonian operator A takes a generating 
function of a conservation law i/A to a symmetry A ( i/) . 

This is the formula for the variational Schouten bracket of two bivectors: 
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= 4i , v/l (A2 ( ^2) ) - 4i ,v/2 (^2 ( n ) ) 
+ ^2 , v/i ( A 1 ( V2 ) ) - 4 2 , V2 1 ( V^l ) ) 

-Ai{B*2{xjfuW2))-A2{Bl{WuW2)), 
where = ^a(va) ~'^^v' '•^^ operators are defined by the equalities: 

iFAi-A*tp=Bi{F,-) onr, 
B*{\l/i,\l/2)=B*'{\lfi,\l/2)\^. 

Here *' denotes that the adjoint operator is computed with respect to the first ar- 
gument. The operators B^ are skew- symmetric and skew-adjoint in each argument. 
Note that if i^" is in evolution form then ( i/^i , Vi) = i^/j ( V^i )• 

Now we are in position to give a definition of a Hamiltonian structure for a gen- 
eral PDE. 

Definition 2. A variational bivectorA is called Hamiltonian if [[A, A]] = 0. 
A Hamiltonian bivectorA gives rise to a Poisson bracket 

{V'i,V'2}a = -E^(v.i)(V^2) + 4*(v^2), (9) 

where \\fi and V2 are conservation laws of <^ and the operator A is defined by the 
relation ^f(A(v/i)) = ^(j^)- 

As in the evolution case, we call an equation bi-Hamiltonian if it possesses two 
Hamiltonian structures Ai and A2 such that [[Ai,A2]] = 0. 

An infinite series of conservation laws y/i, v/2, ■■•is called a. Magri hierarchy 
if for all i we have Ai(y/i) = A2(v/i_|_i). In the standard way one can show that 
{xj/i, v//}ai = {Wii Vj}A2 = for all i and j. 

Now let us return to the question of invariance of the ^* -covering. Suppose the 
equation # under consideration is embedded in two different jet spaces 




We encountered an example of this situation when discussed the KdV equation, with 
7" being jets with coordinates x, t and u, while 7~ being jets with coordinates x, t, u, 
V, and w. Now, we have two linearization operators, l]g and (^g, the former computed 
using the embedding S ^ J" and the latter is obtained using the embedding S — > 7^. 
It is not difficult to show that these two hnearization operators are related by the 
following diagram: 
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— \ 


• 


a 




— > 



(10) 



^2 

where all arrows are '^^'-differential operators on S" satisfying the following relations: 
4j3=j3'4, 4a = a'4, J3a = id+*i4, aj3=id+*24- (H) 



We use the dots • to avoid introducing new notations for the corresponding spaces of 
sections of vector bundles. 

Definition 3. Two '^-differential operators A] and A2 on S are called equivalent if 
there exist ^-differential operators a, ji, a', ji', si, and S2 such that 

AiP=p'A2, A2a = a'Ai, j3 a = id+si^li, aj3=id + S2^2- 

(see [2] and references therein). Thus, we can say that the linearization operators 
and are equivalent. 

The following simple Lemma explains why this notion is really important. 

Lemma 1. ^-differential operators Ai and A2 are equivalent if and only if the Al- 
and A2-coverings are isomorphic as linear coverings. 

So, to prove that ^*-covering is invariant we have to establish that the operators 
i^J and are equivalent. This is impUed by the following result. 



Theorem 3. If two normal operators Ai and A2 are equivalent then A^ is equivalent 
to Al. 

Corollary 1. The equation ^* {S ) does not depend on the embedding S — » 7°°. 



Now, recall that bivectors were defined as conservation laws on J^*{<f?), while 
operators that correspond to them are essentially generating functions of these con- 
servation laws. Thus, the operators depend on using an embedding S J°°. Assume 
that we have two different embeddings as above, so that they give rise to two opera- 
tors A' and that correspond to the same bivector. Here are the formulas that relate 
these two operators: 



A^ = aA^ a'*, 
a1=/3a2/3'*. 



(12) 
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4 Examples 

Let us revise the three examples from the Introduction. 

Example 4 (KdV). We considered two different embeddings of the KdV equation to 
jets: 

Ut-Uxxx- ^UUx = 0, 

Ux-V \ 

Vx — W = 0- 
y — Ut + 6uv J 

Here are all operators of diagram (10): 

= A - Dxxx - 6uDx - 6ux, i% = 



Dx 


-1 







Dx 




A + 6v 


6m 


Dx 





/3 = (1 0) , 
l3'={-Dxx-6u -Dx -l); 




















Dx 


1 






a = \ L>x \ , a' 



Sl = 0, S2 

Formulas (12) relate Hamiltonian operators (1) and (3). 

Remark 5. If we take an operator from (1) for A' and compute via (12) we will 
get an operator from (3) only up to the equivalence. 

Examples (Camassa-Holm equation). The Camassa-Holm equation written in the 
usual form — Utxx — uuxxx — '^UxUxx + 2)UUx = has a bi-Hamiltonian structure: 

Ai=Dx A2 = -Dt -uDx + Ux. 

If we rewrite the equation in the form 

nit + unix + 2uxni = 0, 
m — u + Uxx = 

then the bi-Hamiltonian structure takes the form 
'^l - n _ n3 n ^2 



Dx-D^ OJ ^ \2mDx + mx 



Note that the operators Bi and B2 from Example 2 are entries (up to sign) of the ma- 
trix A\ and Aj. Thus we see that studying bi-Hamiltonian structure of the Camassa- 
Holm equation does not require the use of the (1 — D^)^' "operator". 
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Example 6 ( Kupershmidt deformation ). Let ^ be a bi-Hamiltonian equation given by 
F = and Ai and A2 be the Hamiltonian operators. 

Definition 4. The Kupershmidt deformation SofS has the form 

F+A\{w)=0, Al{w)=0, 

where w = {w^ , . . . , w') are new dependent variables. 

Tlieorem 4. The Kupershmidt deformation S is a bi-Hamiltonian system. 

The proof of this theorem consists of checking that the following two bivectors 
define a bi-Hamiltonian structure: 

(0 ip+Al{w)+A*{w)) ) 

The generalisation of Kupershmidt's theorem from the Introduction is the follow- 
ing. 

Theorem 5. /f i//], \j/2, ... is a Magri hierarchy for S then, under some technical 
assumptions, {xj/i, — y/i+i), i = I, 2, . . . , is a Magri hierarchy for the Kupershmidt 
deformation §. 

Details and proofs of Theorem 4 and 5 can be found in [5]. 
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